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Abstract 

Several lines of reasoning suggest that there might exist a non-sequential 
fourth generation of heavy quarks having very small mixing with light quarks 
and hence exceptionally long lifetime. 
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The accurate measurement in 1989 of the width of the Z boson showed that there 
exist precisely three neutrinos couphng as in the Standard Model (SM) and lighter than half 
the Z mass, ~ 45GeV. 

This led to the conventional wisdom that there are three and only three quark-lepton 
families. The discovery of the top quark in 1995 was thus the final fermion of the 
SM. The only remaining particle is the Higgs boson expected to lie between ~ 60GeV and 
~ eOOGeV. 

This neat picture of just three quark-lepton families and a Higgs as the entire light 
spectrum of matter fields has great appeal. Depending on the mass of the Higgs [Q, it could 
be the entire story up to the Planck or at least the GUT scale. Nevertheless, apart from 
the unsatisfactory aspect that this picture does not explain why there are precisely three 
families there are two principal reasons for suspecting that it is incomplete and that there 
might exist more light particles: 

(1) The strong CP problem is unresolved. Although weak CP violation may be accom- 
modated through the KM mechanism, solution of the strong CP issue is more satisfactorarily 
addressed by spontaneous CP violation in a model with two additional flavors of quark [Q. 

(2) The three couplings of the SM fail to evolve to a common unification point. Until 
recently this was thought |^ to offer support for low-energy supersymmetry, although this 
has been questioned [0] long ago. Very recently one of us [§] has pointed out that a fourth 
family with a Dirac-mass neutrino, because it leads to UV fixed points for the corresponding 
Yukawa couplings, can lead to a satisfactory unification at ~ 3 x lO^^GeV. 

These two quite different reasons lead us to the conclusion that it is quite likely that 
the top is not the last flavor of quark and that there likely exists one further doublet {U, D) 
which is ezt/ier vector-like ((1) above) or chiral ((2) above). 

The mass splitting of this extra doublet is severely constrained by precision electroweak 
data, conveniently parametrized by the S, T variables - the U variable is non-restrictive 
in this context. For the non-chiral case, there is no contribution to S and T at leading order 
unless the doublet has a mass splitting leading to a contribution similar to the first term in 
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Eq(S). 

For the chiral case, there is a contribution to T given by: 
^ lAM^I 1 , |AM|| 1 



(1) 



where Mq is the heavy quark mass, AMg is the mass sphtting in the quark doublet, and 
AM| is the corresponding mass sphtting of the lepton doublet, assuming no Majorana mass 
for Nr. Experimentally T < 0.2 so for any Mq > 200GeV, we deduce that the ratio of the 
U mass to the D mass cannot exceed 1.1. For the chiral case there is also a contribution to 
S given by, for a complete chiral family: 

S=^ = +0.21. (2) 

which is only just compatible with precision data if T < 0.2. 

This leads to our main point: the ratio of masses in the fourth family is 1.1 or less, 
while in the third family we have mt/rrib ~ 30. In the second, rric/ms ~ 10. This suggests 
heuristically that the fourth family is very different, and hence likely to be isolated from 
the first three families by tiny mixing angles. The question then is how tiny this mixing 
might be and how long-lived these heavy quarks can be. For example, let the mixing angle 
between the fourth family U quark and b quark be Vub = x and assume that Vot = Vub- We 
shall particularly be interested in the following two ranges: 10~^ < a; < 10^'^ and x < 10~^. 
The first range is the one considered by models such as the aspon model. For x < 10~^, 
we may assume an almost unbroken 3 + 1 structure under a horizontal family symmetry to 
isolate the fourth family |]10|, similar to the 2 + 1 family structure used in previous models 



TT| which successfully accommodate the top quark mass. 

Without specifying a particular model, we will examine the phenomenology of heavy 
quarks of long lifetime corresponding to small x in the quoted range. We will concentrate, 
in particular, on facilities such as hadron colliders (Tevatron, LHC), where one has the best 
chance of finding these objects. We will first estimate the production cross section for such 
heavy quarks. We then discuss various signatures. 



For the Tevatron, the production cross section can be easily estimated since it is similar 
to the one used for the top quark. The process pp —>■ QQ + X, where Q is a heavy quark, 
can proceed through qq — > QQ and gg — > QQ. At the Tevatron, the qq process dominates 
over the gg one (roughly 90% to 10%). For example, if mq — 180 GeV, the cross section is 
~ 4 pb for ^ = 1.8 TeV and ~ 5.5 pb for = 2 TeV. Also, for ^ = 2 TeV, a heavy 
quark with mass ~ 230 GeV can have a non- negligible cross section of ~ 1.5 pb. Of course, 
the cross section increases tremendously at the LHC, by roughly three orders of magnitude. 
Since the up and down heavy quarks (U and D) are supposed to be fairly degenerate, their 
production cross sections are practically the same: we will have an equal number of UU and 
DD. Their signatures, however, are very different. 

We shall assume that mu > mo- Therefore both U and D can in principle have the 
following decay modes: U — >• or q) + (/"'"t^, 51^2), U ^ q + W, D ^ {t or q) + ^1^2), 
and D ^ {t or q) +W. Here t, q and I denote the top quark, the light quark and the 
light lepton respectively. In those decays, one has, in principle, to distinguish the chiral 
case from the vector-like one. If {U, D) is a chiral doublet, all of the above decays will be 
of the V-A nature. If it is a vector-hke doublet, the process U — > (D or q) + (Z+i/, ^1^2) will 
contain both V-A and V-l-A in the heavy quark current. However, because the light quarks 
are chiral, it turns out that the decay rate is practically the same as if {U, D) were chiral. 
This is because the rate is proportional to (G| -|- Gj^)g1li + GLGRg\l2 and that I2 ~ —h 
and Gl — Gr — gL {h,2 are phase space integrals). The pure V-A case will correspond to 
Gr = and one can see that the two rates are the same. For all other decays involving the 
transition between a heavy and light quark, it will be pure V-A. Therefore, we shall only 
hst formulas related to the pure V-A cases. 

The three-body and two-body decay widths, and are given by 
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FJ =12\Vq^q^\ — — h{mQjmQ^,mw/mQj, 



(3a) 



= ^^^^\VQ,Qj'h{mQjmQ,,mw/mQ,), (3b) 



where I3 and I2 are well-known phase space factors |T2[. Also, Vq^Qj denotes the mixing 
between the two quarks, Qi and Q2- For instance, we may assume that Vjjd ^ ^ and 
\Vub\ ~ \VDt\ = X, where x is to be estimated. We shall start with the decay of the D quark 
first since it will set the range of the mixing parameter x where one can consider at least 
one of the two heavy quarks to be long-lived. We then discuss the characteristic signatures 
for such long-lived quarks. 

The current accessible but unexplored decay length for a long-lived heavy quark to be 
detected is between lOO/im and Im, a range on which we shall focus. (It should be noted 
that decay lengths less than 100/im and greater than Im are accessible as well with the 
latter being explored. Also decay lengths of the order of a few tens of cm might be hard 
to detect.) Since cr ~ 2 x 10^^^ GeV/T{GeV)fim, this would correspond to a width, F, 
between 10^^^ GeV and 10^^^ GeV. (For comparison, the top quark width is approximately 
1.6 GeV.) 

We shall assume that nif) > 175 GeV. 

For the special case of mu ^ 175 GeV, we only have D (c, u) + W. The decay width 
will be Ff ^ 1-6|Vd(c,m)P GeV. This would correspond to a lifetime t ^ 4 x 10'"^^ \Vd{c,u)\~'^ 
s, with cr ^ 1.2 x 10^^°| VD(c,u)|~^At"^- The decay length is given by P'jcr, with /?7 being 
typically of order unity. The present experimental resolution is > 100 fim. For a 175 GeV 
D quark to be detected, the mixing |Vd(c,u)| < 10^^. In the most simple minded scheme for 
the quark mass matrix, one might expect |Vdc| to be at most if {Vntl = x which would 
imply that x < 10^^. This, however, is highly model-dependent. It is best to search for the 
direct decay of D into the top quark. 

When rriD > rrit, D can decay into t via the three and/or two body processes, depending 
on its mass. For rriD between ~ 177 GeV and rrit + my/ ~ 256 GeV, the decay into t is 
exclusively three-body. However, because one can also have D c + W, the total rate 
highly depends on how D mixes with t and c. It is straightforward to compute the width 
for D — > t + {l^u, qiq2) and compare that with D c + W . To do this, let us put \VDt\ = x 
and \Vdc\ = where we shall assume that n >2. 



Let us first assume x < 10~^. The detectability of the D quark in the mass range from 
177 to 190 GeV will depend on what \Vdc\ might be. It goes as follows. Let us first assume 
\Vdc\ = x"^- Let us also distinguish two classes of models: one with x < 10~^ and one with 
10""^ < X < 10"^ (the aspon model is one example of the latter case.) 

In Table 1, we show results for Tf (D ^ T^{D t)/x^ and T^{D ^ c)/\Vdc\'^, as 

functions of mo- We see the following features. I) For the first class of models with very small 
mixing, x < 10"^, it is straigthforward to see that one always has T^{D — > > r^(D — > c). 
Furthermore, one can also see that the decay of D will be detected (within the range of 
1/xm to Im) only when mo > 190 GeV ( corresponding to Fa ~ 2 x 10~^^ for x ~ 10~^). 
The D quark with mass less than 190 GeV will just escape the detector in this scenario. 
As the D quark mass increases beyond 190 GeV, the mixing between D and t will have 
to be smaller if one were to observe the D decay. For example, when niD = 310 GeV, 
Ff (D t)/x^ ~ Ff (L> t)/x^ ~ 3 GeV. This would require x < 10"^. In any case, the 
characteristic signature for this scenario is that D will predominantly decay into the top 
quark! II) For the second scenario with "larger" mixing (such as the aspon model), namely 
10""^ < X < 10~^, there is always an x (e.g. x ~ 10~^) where F^(£) — > c) dominates over 
Ff (£) — > t) for the range uid — 177 — 190 GeV. So in this range, it is possible that D will 
decay dominantly into the c quark! Since F^(Z} — > c) oc x^, this decay is observable. As 
mo > 190 GeV, F|'(I? t) will dominate and the decay is predominantly into t, just hke 
the previous case. However, as one can see from Table 1, this scenario is hard to observe 
(i.e. the decay length is shorter than 1/xm) for itld > 250 GeV (F > 10~^^ GeV). More 
sophisticated silicon pixel detectors will be needed for such a case. 

If \Vdc\ — x^T^fiD ^ t) will always be greater than F^(£) — > c) which means that D will 
decay predominantly into t. The statements made above for the first scenario with x < 10~^ 
still hold in this case. For the second scenario with 10~^ < x < 10~^, the "detectability" of 
the D quark is still possible as long as mo < 250 GeV. However, in contrast with the above 
discussion, the predominant decay mode is into the top quark instead of the charm quark. 

We now turn to the decay of the U quark. Since mu/mD < 1-1 (p parameter constraint), 
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mu—mD < 0.091m(7 < mw unless mu > 890 GeV, a strong coupling scenario no^ considered 
in this paper. The U quark decays into a D via a virtual W, namely U ^ D + ^1^2), 
where / and q are the light leptons and quarks. The U quark can also decay into a light 
quark and a real W, namely U ^ q + W where q — b,s,d and with U ^ b assumed to be 
the dominant transition. Which decay mode of the U is dominant over the other depends 
on how degenerate U is with D and on how small U mixes with the b quark. The results are 
shown in Table 2 where we list as a function of the ratio niD/mu- As for r^(C/ — > 6), the 
estimate is straightforward. We obtain /x^ ~ 1.75 - 4.7 GeV for mu = 180 - 250 GeV. 

For the first scenario with x < 10-^ we obtain {U ^ b) < (1.75 - 4.7) x 10"^° GeV. 
For the second scenario with 10~^ < x < 10~^, we obtain (1.75 - 4.7) x 10~^° <T^{U ^ 
b){GeV) < (1.75 - 4.7) x 10"^. These are to be compared with the results listed in Table 2. 

Unless the U and D quarks are very degenerate i.e. m^j/m;/ > 0.98, the decay mode 
[/ — > D + (Z+i/, 51^2) dominates in the first scenario {x < 10~^). A look at Table 2 reveals 
that the U decay length is much less than Ifim. The signals can be quite characteristic: 
there is a primary decay of the U quark near the beam followed some 100/xm or so later by 
the secondary decay of the D quark. One might see two jets or a charged lepton originating 
from near the colliding region followed by three hadronic jets or one jet plus one charged 
lepton. The reconstruction of the event, if possible, might reveal the decay of a short-lived 
quark (the U quark) into a long-lived quark (the D quark). There would approximately an 
equal number DD produced and hence there might be events where only the decay vertex 
of the D is seen. 

For the second scenario with 10~^ < x < 10~^, we can see that, if niD/mu > 0.97, F^ 
dominates and U will decay principally into b. One would not see the type of events with 
one primary vertex separated by a hundred microns or so from the secondary one as we have 
discussed above. For uid/'itiu < 0.95, the situation is similar to the one encountered in the 
first scenario. 

Thus far, we have assumed mu > mo as in the second and third families. We must not 
exclude the possibility that mo > fnu (as in the first family!). The analysis we have given 
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goes through mutatis mutandis, exchanging t ^ b, etc. A principal difference is that we 
may consider hghter long-hved U quarks {e.g. mu < rrit) than we did D quarks. 

We thank A.T. Goshaw and David Stuart for a discussion of the Tevatron detectors. 
This work was supported in part by the U.S. Department of Energy under Grant No. DE- 
FG05-85ER-40219 (PHF) and under Grant No. DE-A505- 89ER-40518 (PQH). 
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TABLES 

TABLE 1. The values of r§{D t + {lv,qq)/x^, V§{D ^ c + W)/\Vdc\^, and 

(-D — > t + W)/x^ as functions of 





177 180 


190 200 220 


250 


270 


290 


310 


rf/x2(GeF) 


4.3 X 10-11 1.1 X 10-^ 


' 2 X 10-6 2.7 X 10-^ 3.8 x 10" 


-3 2.5 X 10-2 


2.8 


3.1 


2.8 


T^/\VDc\\GeV) 


1.66 1.74 


2.05 2.4 3.2 


4.7 


5.9 


7.3 


8.9 



r^/x'^{GeV) 0.86 1.77 2.86 



TABLE II. The width r^{U ^ D + {Iv, qq)) as a function of mo/mu 



niD/mu 


0.91 


0.93 


0.95 


0.97 


0.98 


n 


5.2 X 10-5 


1.5 X 10-5 


2.8 X 10-6 


2.3 X 10-^ 


3 X 10-^ 
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